The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

Generalized Jacobi spectral Galerkin method for fractional pantograph differential equation

This work is concerned with the extension of the Jacobi spectral Galerkin method to a class of nonlinear fractional pantograph differential equations. First, the fractional differential equation is converted to a nonlinear Volterra integral equation with weakly singular kernel. Second, we analyze the existence and uniqueness of solutions for the obtained integral equation. Then, the Galerkin method is used for solving the equivalent integral equation. The error estimates for the proposed method are also investigated. Finally, illustrative examples are presented to confirm our theoretical analysis.     

High-frequency lengthwise diffraction by the curve separating soft and hard part of the surface

The paper examines the model problem of high-frequency diffraction by a convex surface consisting of two parts. One is soft, the other is hard. The incident wave falls at a small angle to the line which separates soft and hard parts of the surface. The change in the boundary condition provokes the field in the Fock zone to have a rapid transverse variation. This causes a special boundary-layer to be formed. The boundary value problem for the three dimensional parabolic equation is reduced to the Riemann problem solved by the factorization in the form of infinite products containing the zeros of the Airy function and zeros of its derivative. the results of this factorization appear under the sign of double Fourier integral in the representation of the field. Both numerical and asymptotic analysis of this representation is carried out and illustrates the effects of high-frequency diffraction caused by the line of the boundary condition discontinuity.     

The Bergman kernel for the Vekua equation

The paper extends some well-known results for analytic functions onto solutions of the Vekua equation $\partial \overline{{}_z}W=aW+b\overline{W}$ regarding the existence and construction of the Bergman kernel and of the corresponding Bergman projection operator.     

POD enriched boundary models and their optimal stabilisation

Modes obtained using the Proper Orthogonal Decomposition are used as boundary enrichment functions within a variational multiscale method for the stochastically forced Burgers equation. Initially, large increases in accuracy are obtained using the enrichment functions without stabilisation terms. Then, optimal coefficients for the stabilisation parameter τ of the unresolved scale model are calibrated using a goal‐oriented model‐constrained optimisation technique, resulting in further improvements. As both the determination of the enrichment functions and the optimisation of the coefficients requires high‐accuracy reference data, a scaling procedure is introduced to allow their use over range of conditions. Numerical experiments confirm that the scaling procedure is effective. Copyright © 2014 John Wiley & Sons, Ltd.